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The doped Mott insulator in one dimension has been studied based on the phase Hamiltonian 

^ i_ 

, with the Umklapp scattering process, in which the charge degree of freedom is described by 

the quantum sine-Gordon model. The well-known equivalence between the quantum sine- 
Gordon model and the massive Thirring model for the spinless fermion makes it clear that the 
Mott- Hubbard gap originates from the Umklapp scattering process as was indicated by Emery 
and Giamarchi. Compressibility density-density correlation function, frequency dependence of 

■ optical conductivity and Drude weight have been calculated in the presence of the impurity 
, scattering treated in the self-consistent Born approximation. It is seen that there exists a 

qq , crossover behavior in the spectral weight of charge excitations: the acoustic mode is dominant 

in small wave number region while the optical excitations across the Mott-Hubbard gap lie in 
' large wave number region and that this crossover wave number is reduced as the Mott transition 

' is approached. 
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§1. Introduction 

The Mott transition had so far been generally studied as metal-insulator transition which is caused 



by changing the strength of interaction (or the band width)Hro ) at a particular carrier density of 
half-filling. Since the discovery of high- T c material, however, the behavior as the band filling (or 
doping rate) is varied in the vicinity of Mott insulator is the central issue, i.e., doping into the 
Mott insulator (doped Mott insulator). The high-T c superconductivities can be considered as one 
of ground states of this doped Mott insulator. In the metallic state of this doped Mott insulator, 
anomalous behaviors are observed in transport and magnetic properties and then this state is called 
anomalous metallic state. The origin of those anomalies are considered to arise from characteristic 
nature of electronic states near the Mott transition. 

For theoretical studies on this anomalous metallic phase, the Hubbard model has been frequently 
used. Although this model is considered to be the simplest Hamiltonian containing the essence of 
Mott transition, little is known about electronic properties near the transition as the doping rate 
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is varied. In one dimension, the Hubbard model has been solved exactly by use of Bethe Ansatz 
H and anomalous behaviors of various quantities, compressibility, Drude weight etc., have been 
disclosed .i'i) Further development has been achieved recently by the conformal field theory and 
asymptotic form of correlation functions have been determined.!' These correlation functions 
behave as the power law which is characteristic to the Tomonaga-Luttinger liquid! 0) where there 
exist excitations of acoustic modes both for charge and spin degrees of freedom in the limit of long 
wavelength. However, the detailed nature of the excitation spectrum has not been explored with 
this framework. 

On the other hand, these properties can be explored by the bosonization.0>lll>lll) The system 
far away from half-filling is in the Tomonaga-Luttinger regime. As for a Mott insulator, i.e., 
at half-filling, Emery has shown that the Mott-Hubbard gap is seen to result from the Umklapp 
scattering process which is neglected in the Tomonaga-Luttinger regimeB© The Mott transition 
caused by carrier doping, i.e., the doped Mott insulator, has been studied by Giamarchi.0 In this 
study, optical conductivity was calculated within the second order of Umklapp scattering process 
by use of memory- function approximation and renormalization group method.ll! 00) The Drude 
weight and its critical behavior toward the Mott transition has been investigated by Giamarchi@) 
and Emery.© These studies on the optical conductivity are concerned with charge excitation at 
q = 0. On the other hand, charge excitation for finite q and oj near the Mott transition has been 
disclosed in our previous work@ where we have found interesting features of spectral weight of 
charge excitations; i.e., the crossover behavior where the acoustic mode, which is the property of 
Tomonaga-Luttinger liquid, has a large spectral weight in the region of the small wave numbers 
while the optical mode, which originates from the inter-band excitation across the Mott-Hubbard 
gap, is weighted in the region of large wave numbers. The characteristic wave number of this 
crossover is reduced to zero as the metallic state merges to the insulating state. 

In this present work, we will study in detail effects of disorder on charge excitations in the doped 
Mott insulator. To investigate such effects, the phase Hamiltonian0'llIlllI|) has been employed, 
where the charge degree of freedom is described by the quantum sine-Gordon model (QSG model) 
with misfit parameter representing the proximity to the half-filling. §13) The phase Hamiltonian has 
been originally derived by the bosonization to describe the low energy excitations for correlated 
systems in the limit of weak coupling!!) and has been considered to be qualitatively applicable 
to the regime of strong coupling by choosing properly the constants of the Hamiltonian.iHil) In 
the classical treatment of this model, it is known that discommensuration lattice (soliton lattice), 
which will be identified with particles corresponding to the charge degree of freedom, is created 
as the misfit parameter is varied and that there exists a transition between the commensurate and 
incommensurate phase as clarified by McMillan© 

We will study quantum mechanically the spectra of charge excitation by mapping QSG model onto 
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the massive Thirring model (MT model).il'@'i3ilii) This useful mapping was originally applied 
by Luther and Emery to the study of spin dynamics in ID electron systems with backward scattering 
process!]'© and later to the Mott transition@>00@@'S) as mentioned above. However, the 
detailed study of charge excitation including disorder has not been carried out especially for the case 
where the Mott transition is approached by the carrier doping. By making use of the mapping and 
including disorder in terms of spinless Fermion, we will study the compressibility, density-density 
correlation function and optical conductivity in the presence of impurity scattering. 

The structure of this paper is as follows. To make this paper complete, we summarize the results 
of some previous works by Emery, Giamarchi and ourselves in §2 and §3. In §2, we introduce 
the massive Thirring model to describe the charge degree of freedom. This model is derived from 
the phase Hamiltonian by using the bosonization. In §3, the charge excitation spectrum and 
optical conductivity are calculated in a clean system. In §4, we introduce an impurity scattering to 
investigate these quantities in disordered systems. §5 is for conclusion and discussion. We confine 
ourselves to absolute zero temperature in this paper. 

§2. Hamiltonian and One Particle Properties 

2. 1 Phase Hamiltonian and Classical Picture of Excitations 

In one dimensional electron system, the low energy excitation of charge and spin degrees of free- 
dom are independent of each another (spin-charge separation). Therefore, the phase Hamiltonian 
comprises two parts as follows; 

H p = Jdx [A p {V9{x)f + B p cos{29{x) - q x) + C p P 2 ] , (2.1) 

H a = Jdx [A a (Vcj)(x)) 2 + B a cos(20(x)) + C a M 2 ] , (2.2) 

where 

[9(x),P(y)] = iS(x-y), (2.3) 

[<Kx),M(y)]=iS(x-y). (2.4) 

9{x) and 4>{x) describe the fluctuation of charge and spin degree of freedom, respectively. The 
parameters A p , B p , C p , A a , B a and C a depend on the coupling constants and the filling of band. 
The charge excitations, which are our main interest in this paper, are described by H p , eq.( |2.l| ). 
where the .Bp-term is due to Umklapp scattering with q$ = G — 4/co, the misfit parameter. Here, G 
is the reciprocal lattice vector and ko = h/vfo, with /x and vfq being the chemical potential and 
the Fermi velocity of the noninteracting electrons, respectively. In the case of the noninteracting 
system, ko is the Fermi momentum. This misfit parameter controls the doping rate, i.e., the number 
of hole. 

In order to understand the physical implication of QSG model, eq. (|2.1| ), as a function of the 
misfit parameter, we first analyze the model semi-classically by assuming that the field operator, 
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9(x), consists of classical solution, 9 c (x), and fluctuation around this solution, 9(x,t), as 9{x) = 
6 c {x) + 9(x,t). Here 9 c (x) and 9(x,t) satisfy the following equations, respectively, 

A p (V9 c (x)) 2 + B p cos(20 c (x)) = D, (2.5) 

VV(x) + Um 2 + Cj 2 ) - 2m 2 sn 2 (mx/k; k)] cp(x) = 0, (2.6) 
where we assume 9(x,t) = (p(x)e lult and 

2 2B P 

m - 



k 



2 



2m 1 



m 2 + 2D/Ap' 

2 



In eq.(2.6), sn 2 (mx/k;k) is the Jacobian elliptic function with the parameter /c. An integral con- 
stant, D, is to be evaluated by the solution, 9 c (x), which minimizes the energy per unit length, e, 
given by, 

r 2irl 



l - ^ dxAp(V9 c (x) + q ) 2 + Bp cos(20 c (x)), (2i 
vrt Jo 



2ttZ 

where I is defined as, 2-7r/ = (2k/m)K(k), and identified to a mean distance between discommen- 
suration(DC)s. K(k) is the first kind complete elliptic integral with parameter k. In the classical 
solution, 9 c {x) = ir/2 — sn(mx/k;k), there exist a transition between the commensurate(C) and 
incommensurate(IC) phase as is shown in Fig.l for m = 1 and some choices of k. Since the 
fluctuation of number density of charge, n(x), in the long wave length limit (q <S 2ko) is , 

n(x) —, (2.9) 

7T 

the DCs of 9 c {x) are identified as particles.© Thus the commensurate (incommensurate) phase 
can be considered to be an undoped (doped) state. In the IC phase near the Mott transition, the 
DC lattice is created and its excitation spectrum, which is the eigenvalue of eq, (|2.6j ) and is shown 
in Fig. 2 for m = 1 and k 2 = 0. is considered to be vibrations of DC lattice; the acoustic 

mode of the DC lattice in the first Brillouin zone and the optical mode in the extended zone 
In this classical treatment, there exists a critical wave number, ir/l, which separates the region 
of the acoustic and optical mode. This ir/l characterizes how the metallic state merges into the 
Mott insulator. At the C-IC transition, the acoustic mode disappears and only the optical mode 
survives. These characteristic features are held even in fully quantum mechanical treatment, as 
we will see. Although general features of excitation spectrum can be understood by the classical 
picture, the details of this spectrum, the compressibility, and the optical conductivity should be 
disclosed in the quantum treatment. 
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2.2 Massive Thirring model 

To solve eq. (|2.1| ) quantum mechanically, we describe the charge degree of freedom in terms of 
spinless Fermion0@lil'lll'0 instead of Boson. Eq.(|2.1|) is mapped onto the following Hamiltonian 
(See Appendix A), 

H p , F = v c J dx(^(x)(-idT 3 )t>(x)) + ^^(x)^(x) 
+ V J dx^ (x)ti^(x) 

+ dx[(^(x)V(x)f - (*t( z ) ri *(- r ))2] j (2.10) 

where 




W = irA p - 

and Tj, (j'=0,l,2,3) are Pauli matrices. The second term, which originates from the Umklapp 
scattering process, creates the charge gap (Mott-Hubbard gap). This field operator, \&(x), is 
identified as a charge particle because *$>(x) creates a DC in 9{x) as follows!!'©, 

[0(x),ip}(y)} = ^sgn(x-y)rp](y) (j = 1,2). (2.11) 

As mentioned in the preceding section (see eq.(|2.9|)), this DC corresponds to a charge particle. In 
the following discussion, we will set W = and then eq.(2.10) represents a system of free Fermions. 
This particular choice, W = 0, is a good approximation near the Mott transition as has been noted 
by Giamarchi.lll) The reason is as follows; if W = we see 

Ap = ±kk; 

C p = TTVpKp, 

therefore, 

K P = \, (2.12) 

where v p is the velocity of charge excitation and K p is the critical exponent for the density cor- 
relation function. On the other hand, it is known that both quantities, v p and K p , are generally 
determined by the degree of the filling of band and the coupling constant. Especially, K p approaches 
1/2 near the half-filling for any values of coupling const ant .i'lBH* It is noted that K p also becomes 
to 1/2 toward the infinite repulsive interaction for any filling as noted by Emery© 
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2.3 Energy Spectra and Compressibility 

With W = 0, we can solve eq.(2.10) exactly and get a simple physical picture about the charge 
degree of freedom near the Mott transition.^ MEM) First, the energy dispersion is 



E-(k) = ^-y/vlk* + V*, (2.13) 



E+(k) = ^ + yjv*k 2 + V 2 . (2.14) 

The band structure relative to the chemical potential is shown schematically in Fig. 3. The actual 
carrier density (doped hole density) away from half-filling, 5 , is determined by the range of k which 
satisfies E_(k) > 0. For q > q c = 2V/v c , 5 is determined by, 



2lT y v % 2 - (2-15) 

and 5 = for qo < q c . This relation determines the carrier density as a function of the chemical 
potential, go- Consequently we immediately see that the compressibility, k = dn/d/j,, diverges 
inversely proportional to the doping rateES as found in the study by Usuki, Kawakami and Okiji 
based on the Bethe Ansat2&, because of the particular energy dependence of the density of states 
at the band edge in one dimension. 

§3. Charge Excitation and Conductivity in the Clean System 

We will first calculate the density-density correlation function and next the conductivity and 
Drude weight by use of the Nambu formalism (See Appendix B). 

3.1 Charge Excitation Spectrum 

The charge excitation spectrum is given by the imaginary part of the retarded density-density 
correlation function, ImN' R (g,a;). Since we are interested in the spatially slowly varying part (q <C 
2&o), the charge density is given by n(x,t) = e^ (x,t)^(x,t). After straightforward calculation, 
we obtain the following result, 

JmN R (q,u) = — J —J -(/( e )-/(e + ^))Tr[T ImG H (k + q, e + a;)roImG H (k,e)] 

2e 2 fdk r/ V 2 . e . 

= -T- / — [{VkVk+q + u k u k+q + )d(LO + E k+q - E k ) 

n J It: ZE k E k+q 

V 2 

+ (v k u k+q + u k v k+q - )8{u) - E k+q - E k )}, (3.1) 



where tq is the unit matrix and 



1. k . 
Vk = x(l - -=r), 



E k = Jv 2 k 2 + V 2 . (3.2) 
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In the last equation, the k-integration is carried out for v c qo/2 < Ef. < v c qo/2 + u, /(e) is the Fermi 
distribution function and G(k,ie) is Fourier transform of Green function defined as follows, 



G(x,t) = -{T t 



V>i(x,t)4(0,0) V 2 (x,r)4(0,0) 
V>i(x,t)4(0,0) ^ 2 (x,r)4(0,0) 



The detailed explanation of the Green function is given in Appendix B. From eq.(3.1), we see that 
the charge excitation exists in the following regions; 
1. Acoustic Excitation: 

(v c qf + AV 2 - u? > 0, 
(v c q) 2 - uj 2 > 0, 



v c qo -io < {v c q)\ -, ^ 5 < v c q +lo. 



2. Optical Excitation: 



(v c q) 2 + 4V 2 - to 2 < 0, 
(v c q) 2 - u 2 < 0, 



(v c q) 2 - UJ 2 

or 



r ^ \ u cV t<±v —uj 
{v c q - to < v c q\j — — < v c q + to 



{v c qf+W^-^ 

-{v c q -u)> v c q x -5 = > -{v c q + u))\ . 

[VcQ) — w 



These regions are plotted for several choices of doping rate in Fig. 4 and actual magnitude of the 
spectral weight, ImN^(g, uj), in these regions are shown in Fig. 5.© The existence of the acoustic 
mode of excitation spectrum has also been obtained from a quantum Monte Carlo simulation 
using the maximum entropy method for the ID Hubbard model away from half-filling.il) In the 
metallic state (far away from half-filling) , the acoustic mode, which corresponds to the charge 
excitation in the Tomonaga-Luttinger liquid, dominates spectral weight especially around the small 
wave numbers. Thus it is reasonable to consider that the inter-band excitation across the Mott- 
Hubbard gap does not contribute to the density-density correlation function in the limit of long 
wave length. Near the Mott transition, however, such region where the acoustic mode is dominant 
becomes small and the optical mode is dominant elsewhere. This implies that, in the real space, 
the density-density correlation function has a characteristic length scale for the crossover behavior. 
The similar crossover behavior is numerically found in the spin-spin correlation functionJH' As 
the Mott insulator is approached, the characteristic length scale becomes large, above which the 
correlation function follows its asymptotic form derived by the conformal field theory. In the 
limit of Mott insulator, this length scale diverges and the conformal field theory breaks down. 
Correspondingly, the spectral weight shifts from the acoustic mode to the optical mode and the 
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metallic phase crossover to the insulator phase. The Mott transition is characterized by this 
crossover,© which has not been so far found out by any theories. The location in g-space of 
this crossover about the dominant spectral weight corresponds to that of the excitation spectrum 
in the classical approximation as shown in Fig. 2. 

3.2 Conductivity and Drude Weight 

Next the conductivity@0>lll) and Drude weightEBIll'il) will be calculated by the imaginary 
part of retarded current-current correlation function, ImK^(w), which is evaluated as follows. By 
noting that the current density is j(x,t) = ev c ^'(x, i)r3^(x, t) because of the continuity equation, 
we obtain 

9p^i^ r c\k P r\f R R 

lim ImK^, u) = ^ / ^ / —(/(e) - f(e + w)) lim TV [ImG> + q,e + u)r 3 ImG H (k, e)r 3 ) 

q^O n J ZlT J 2lT q— >0 



q 



2e 2 v 2 f dk V 2 

lim[(v k v k+q + u k Uk+q - )8(u) + E k+q - E k ) 



h J 2nq^o iK ^ 2E k E k+q 
V 2 

+ (v k u k+q + u k v k+q + )S(u - E k+q - E k )\. (3.3) 

^^k-^k+q 

In the last equation, the k-integration is carried out for v c q^j2 < E k < v c qo/2 + to. 

Near the Mott transition, the Drude weight, D c , which is defined as, Rec"(u;) = D c 5(uj), is calculated 

as follows (see also Appendix C), 

„ , , „ K r (lu) -2v c e 2 /h , ^ 

Rea(uj) = Re — ^— i— , (3.4 

ILO — 7] 

therefore, 

e 2 v c 5 



'c 



n go/2' 

e 2 v r 5 



n v < 3 ' 5 ' 

This fact, that D c is proportional to the concentration of the doped holes near the Mott transition, 
is consistent with the exact result obtained by Shastry and Sutherland by use of the Bethe Ansatz 
I!* and has also been noted by Giamarchill!) and Emery.ll!' It is interesting to notice that the 
doping dependence of D c is different from that of Tomonaga-Luttinger liquid. In the Tomonaga- 
Luttinger regime, the D c is determined as, 2e 2 v c /h, by the electron (particle) number density which 
is expressed through This result is obtained by the Kubo formula by using the anomalous 

commutation relation which is inherent to the one dimensional systems with linear dispersion.!!^ 
The different dependence on the carrier concentration of D c between in the Tomonaga-Luttinger 
regime and near the Mott insulator results from the roles played by the Umklapp scattering process. 
Thus we can understand that the Umklapp scattering process plays crucial roles near the Mott 
transition. 
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§4. Impurity Effect 

So far we confined ourselves to the clean system. We will in the following consider the effects 
of disorder represented by the impurity scattering. We first explain how to introduce and treat 
impurity scattering, and then we study the effects of impurity scattering on the density of states, 
which determines the compressibility and the conductivity. We will employ the self-consistent Born 
approximation and ignore the possible consequence of the Anderson localization or the bound states 
in the Mott-Hubbard gap in order to understand overall features in the present study. 

4-1 Impurity Hamiltonian and Self-Energy 
The impurity potential, u(x), couples to the charge density, p(x), as follows, 

H imp = ^2 J dxu(x - Ri)p(x) (4.1) 

i 

where p(x) is expressed as follows, 

p(x) = 1p\{x)lpl{x) + 1pl(x)lp2{x) 

+ e-' lXx 4 (x)^ 2 (x) + e iXx 4 (x)V>i (x) . (4.2) 

The oscillating part of p(x), which is characterized by the constant, A, originates in the density 
fluctuation across the Fermi surface of spinless fermion. We suppose the impurity potential to be 
the ^-function, u{x) = u5(x). In this impurity scattering Hamiltonian, we ignored the coupling 
between the charge and spin excitation through the 2/^F-oscillating part of charge density© 

In the second order self-consistent Born approximation, the self-energy, which is diagrammatically 
shown in the Fig. 6, is given as follows; 

a(z) At(z) 
A(z) a(z) 

li(k,z) + g 2 (k,z) ?Hk,z) 

f(k,z) gi{k,z) +g 2 {k,z) 



t(z) 




(4.3) 



where 



n-iU 2 
a{z) = -l- 



v c ^J~ z 2 _ y2 

A(*)=AT(* =-i-J — 4.4 

2v c ^Ji2 _ y2 

which are easily obtained after k-integration. For details of z, V, see Appendix B. The self-consistent 
equations eq.(4.3) and eq.(4.4) are reduced to the following equations, 

UiU 2 z 

z = z — a(z) = z + i- 



V c ^J~ z 2 _ y2 

V = V + A(z) = V-\ n%u2 r J ^— . (4.5) 



9 



Eq.(4.5) is rewritten as 



X 4 - 2y X 3 + (y 2 + C 2 " 1)X 2 + 2y X -y 2 = 



(4.6) 



where (, = 3niu 2 /2v c V , x = z/V and y = z/V. 



4-2 Density of States and Compressibility 



By using the solution of eq^^), the density of states, T>(e), is calculated as follows, 



1 r r\k 

V{e) = — / 2-ImTr[G fi (M], 

7T J ZTT 




(4.7) 



T>{e) is plotted in Fig. 7 for several choices of dimensionless parameter, £, representing the strength 
of impurity scattering. It should be noted that, in the presence of impurity scattering, the Mott- 
Hubbard gap, which is considered to be the energy gap in the density of states of spinless Fermion, 
exists for < £ < 1. In the following, we will use the term 'weak disordered system' to refer to the 
system for < £ < 1. On the other hand, for £ > 1, the gap disappears but the density of states 
has a cusp at e = 0. 

In addition to the change of Mott-Hubbard gap, the V{e) in the disordered systems does not 
diverge at the band edge in contrast to the T>(e) in the clean system. This qualitative change of 
density of states directly affects the compressibility, k = dn/dfj,; In the clean system k diverges, but 
in the weak disordered system, it smoothly vanishes as <5 2//3 as is shown in Fig. 8. Though the value 
of the critical exponent 2/3 should not be taken literally in view of the approximation employed, 
this drastic qualitative difference between the clean and disordered systems is to be noted. For 
£ > 1, k approaches a finite constant. 

4-3 Charge Excitation Spectrum in the Disordered System 

In the calculation of the density-density correlation function, N(q,uj), we need the vertex correc- 
tion in addition to the self-energy correction. The vertex corrections to the density vertex, T^(q, lo), 
is the solution of following Dyson equation, 



T N (q, iuj) = t + — J — \G(k + q,z + iu;)r G(/c, z) 
+ r + G(k + q,z + ioj)T G(k, z)t~ 
+ r~G(k + q,z + iio)T G(k, z)t + 
+ G(k + q,z + iuj)T N (q, iu)G(k, z) 
+ T + G(k + q,z + iu)F N (q, kj)G(k, z)t~ 
+ T~G(k + q, z + iw)rAr(g, iu>)G(k, z)t + , 




(4.8) 
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where = [t\ ± ir 2 )/2 and tq is the unit matrix. N(q,iu>) is calculated by using the solution of 
the above equation as follows, 

N(q,iw) = -(-T)^2 J — Tr [r G(A; + q, z + iu)T N (q, kj)G(k,z)\ . (4.9) 

ImN(q, uj) are plotted in Fig.9 for several values of 5, ( and q, by solving the above equation 
numerically. 

In the weak disordered systems, there also exist the acoustic excitation and the optical excitation as 
in the clean system, though these excitations are damped. Hence, the crossover feature of spectral 
weight in q and cu-space survives as far as the impurity scattering is not so strong as seen in Fig.9 
for ( = 0.1 and 5 = 0.3. The doping dependence of spectral weight is shown in Fig.9 for £ = 0.1 and 
q = 2, where it is seen that the dominant excitation smoothly changes from the acoustic (optical) 
excitation to the optical (acoustic) excitation as the doping rate is varied. Therefore the crossover 
behavior in the vicinity of the Mott transition is visible not only in the clean system but is present 
also in the weak disordered systems as well. 

4-4 Optical Conductivity and Drude Weight in the Disordered System 

In the calculation of the conductivity, we need the vertex correction for the current vertex, 
Ta"(0, to), again similar to that of density-density correlation function; i.e., insert T3 between Green 
functions for the first, second and third term in eq.(4.8) and take the limit of q — > 0. 

2C f Ak r - 

r K (iu>) = t 3 + ^ / — |G(fc, z + i^)r 3 G(fc, z) 
+ r + G(k, z + iw)r 3 G(/c, z)t~ 
+ T~G(k, z + iuj)T 3 G(k, z)t + 
+ G(k + q,z + iuj)t K (iuj)G(k, z) 
+ r + G(k + q,z + iw)r K (ku)G(k, z)t~ 

+ r~G(k, z + iu)f K (kj)G(k, z)t + ] . (4.10) 
The current-current correlation function, K(kj), with vertex correction is calculated as follows, 

K(ic) = ^(- T)E / ^Tr [t3G(*, z + iu)t K ^)G(k, z)] . (4.11) 

The optical conductivity is obtained by solving the above equation numerically. In Fig. 10, the 
w-dependence of optical conductivity in the case of 5=0.01 is shown for several choices of £=0.0, 
0.01, 0.1, 0.3. There exist optical absorption across the Mott-Hubbard gap as well as the Drude 
tail. Because of the impurity scattering, both spectra become broad and difficult to distinguish 
from each another if the impurity scattering becomes strong. An example of the doping dependence 
of the optical conductivity in the weak disordered systems is plotted for ( = 0.1 and some choices 
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of 5 in Fig. 11. In this figure, the Drude weight seems to behave similarly to that of clean system 
as a function of doping rate, i.e., smoothly approaches to zero. It becomes clear by comparing the 
integrated weight of both, Drude and optical, that the missing weight of Drude peak shifts to the 
optical absorption and the metallic state merges to the insulator phase. These general features of 
optical conductivity is plausible but it has not been explicitly demonstrated so far. The crossover 
feature near the Mott transition appears as the smooth shift of weight in the optical conductivity. 

§5. Conclusion and Discussion 

In the present study, we investigated the excitation spectrum of density-density correlation func- 
tion and the optical conductivity near the Mott transition in one-dimensional electron system in 
the presence of disorder. Those quantities were calculated based on the massive Thirring model for 
the spinless Fermion. In the Mott insulator, i.e., at half-filling, the Mott-Hubbard gap was seen to 
originate from the Umklapp scattering process as noted by Emery. DDE Even in the doped Mott 
insulator, i.e., away from half-filling, this process played the important role as indicated by Gia- 
marchi.0) First, we calculated the static limit of the density-density correlation function and the 
conductivity, i.e., compressibilityDD 1 and Drude weig In the clean system, the doping 

dependence of both quantities were same as the results deduced by the exact solution.!* Although 
the Drude weight did not qualitatively change its behavior also in the weak disordered systems, 
the compressibility even in the weak disordered systems behaved completely different way from 
that in the clean system. Next, we studied the dynamical properties. In both the clean and weak 
disordered systems, the crossover behavior existed in the spectra of density-density correlation 
function; in the region of small wave number, the acoustic mode was dominant while in the large 
wave number region, the optical mode was weighted. By decreasing the doping rate, the region, 
in which the acoustic mode was dominant, became small and its spectral weight shifted to that 
of optical mode. We would like to emphasize that the Mott transition was characterized by the 
crossover feature rather than drastic phase transition. The smooth shift of spectral weight from the 
Drude weight to the optical absorption is present in the optical conductivity in both the clean and 
weak disordered systems. If the impurity scattering becomes strong and the effect of localization 
becomes important,© the optical conductivity would vanish as approaching uj = and will show 
an characteristic features of the pinning frequencyS 
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Appendix A: Mapping from phase Hamiltonian to massive Thirring model 

The equivalence between the two Hamiltonians eq.(2.1) and (2.10) is proved by the bosonization 
method. To begin with, we define the boson operators, 6j(x) (j = 1,2), in terms of the spinless 
Fermions ipj(x) (j = 1,2) as follows, 



(x) = ~(pi(x) + e 2 (x)), 



2tt 1 
9 j = 1 2l j^exp(-iqx - -a\q\)pj(q), (j = 1,2) 



(A-l) 
(A-2) 



k 

Pj(-q) = Yj 4,k c J,k+ q , (q > o), 



where 



[pi(-q),pi(q')] = [p 2 (q),p 2 (-q')] 



Lq 
2vr 



$q,q'- 



The spinless Fermions correspond to those Boson operators as follows, 

1 



V2 



ira 



: exp(±i0»(:c)), (+ for i = 1, - for i = 2). 



(A-3) 



(A-4) 



(A-5) 



The Hamiltonian eq.(2.1) is described in terms of the density operators as follows, 

/ dx(V(9(x)) 2 = vr / dq [pi(q)pi{-q) + p2{~q)p2{q) + Pi(q)p2(-q) + Pl(-q)p2(q)] , 

J Jq>0 

f dx(P{x)) 2 = - [ dq [px{q)px{-q) + p2(-q)P2(q) - Pi(q)P2{-q) - Pi{-q)P2{q)\ 

J VI" Jq>0 



thus, 



A P J dx{V9{x)f + C p J dx(P(x)) 2 =(^A p + ^j J >o dq\p 1 (q)p 1 (-q)+p 2 (-q)p 2 (q)} 



q>0 

| KA p - j dq [pi(q) P 2(-q) + Pl(-q)p2(q)] (A-6) 

7T / Jq>0 



Next, the first and third term of eq.(2.10) are also written in terms of the density operators as 
follows, 



(A-7) 

(A-8a) 
(A-8b) 



H l% = v c dq [pi{q)pi{-q) + pa(-g)A»a(g)] , 
H f ] D = W I dq [pi(q)p 2 (-q) + Pl(-q) P 2(q)] 



lq>0 

which are proved by the following equations, e.g., 



H{ P %ipiAq) 



v c qpi,a{q), 
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for the first term of eq.(2.10). By comparing eq.(A-6) and eq.(A-7), we see the following relations, 

v c = irA p + C p -, 
W = ttA p -C p -, 

B p term is easily derived by inserting eq.( |A-5| ) into the second term of eq.(2.10). 

Appendix B: Definition of Green Function 

In the Nambu formalism for the charge degree of freedom near the Mott transition, the Green 
function, G(x,t), is defined as 2 x 2 matrix, 

G(x,t) = -(T t \ , . 

V ^i(x,r)^ 2 T (0,0) V 2 (x,r)4(0,0) 

(B-l) 

where T T means the chronological order in the imaginary time. In the clean system, after Fourier 
transformation for the system without the interaction, W = 0, 



G°(k,z) 



f°(k,z) g° 2 (k,z) 

z — v c k —V 
—V z + v c k 



(B-2) 



where z = ie — v c qo/2. 

In the disordered systems, full Green function, G(k,z), including the self-energy correction is 
calculated as follows, 



-i 

(B-3) 



G(k,z) = f&^z^-tizfj 1 , 

z -v c k-a(z) - V - At (z) 
-V - A(z) z + v c k-a(z) 

We take the full Green function as follows to determine the self-energy (or Green function) 

gi(fc,z) ft(M) 

f{k,z) g2(k,z) 



G(k,z) 



z — v r k —V^ 



-V z + v c k 



(B-4) 

We can derive the self-consistent equation eq.(4.5) by substituting eq.(B-4) to eq.(4.3). 
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Appendix C: Calculation of Drude Weight near the Mott Transition 

An uniform electric field, E x (t), along the one dimensional system has such a relation, E x (t) = 
<9 x Ao(x,t)— <9 t Ai(x,t), to the vector potential, A(x,t) = (Ao(x,t), A\{x,t)). To make our discussion 
simple, we introduce external field as follows, A(x,t) = (xE x (t),0). The perturbation term by the 
external field, Ti. ex t, is introduced as ordinary way, 

H cxt = - J dxn{x)A {x,t). (C-l) 

On these grounds, the conductivity, a(uj), is calculated as follows, 



<[J(0,0),n(0,y)])y 

/oo 
dy([J(0, S ),J(0,y)]). (C-2) 
-oo 



h 


roo roo 

/ ds / dy 

JO J-oo 




roo _l 


= h 


/ dy-. ( 


/_oo 1L0 - 7] 




roo J\oj-ri)s 


+ I 


/ ds 




JO \UJ - 7] 



where the translational invariance in space coordinate and the continuity equation are noted. 
In eq.(C-2), the first term in the last equation is estimated as, 

i roo _i 2(= 2 7) —1 

-/ dy-. ([J( ,0),n(0,y)])y = — ^, , (C-3) 

fi J-oo \UJ — T] h VjJ — T] 

by using the anomalous commutation relation, 



w, 



7T 

As a result, the real part of conductivity is calculated as follows, 

K R (u) - 2v c e 2 /h 



[J(x),n(v)] = —9x6(x-y). (C-4) 



Rea(uj) = Re- 



1UJ — 7] 

R 2v c e 2 \ ImK R (w) 



= <^ ReK K (cj) - ^— (-7r)<5(w) + ^. (C-5) 

I h J LO 

Without the Umklapp scattering process in the last equation, the terms corresponding to the K R {u) 
is zero and the Drude weight takes the constant value, 2e 2 v c /h. 
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FigureCaption 

Fig.l The classical solution is plotted for m=l and k 2 =0.5, 0.9, 0.99, 0.999. 

Fig.2 Energy spectrum of the fluctuation around the classical solution is plotted for m=l and 
fc 2 =0.5. 

Fig.3 The band structure of the spinless Fermion (the charge degree of freedom for finite S. 
Fig.4 The region of (q,u) with a finite lmN(q,u>) for several choices of (a) 5 = 0.7, (b) 5 = 0.3, (c) 

5 = 0.01, (d) 5 = 0.0 per one site. 

Fig.5 The spectral density of the charge excitations, lmN(q,uj), in the plane of q and u for (a) 

6 = 0.3, (b) S = 0.1, (c) 6 = 0.0 per one site. 

Fig.6 The Feynman diagram for the self-energy. The thick line means the full Green function and 

the broken line represent the impurity scattering. 

Fig.7 The density of states is plotted for C=0, 0.01, 0.1, 0.3, 0.5, 1.2. 

Fig.8 The compressibility is plotted for C=0,0.01,0.1,0.3,0.5. 

Fig.9 The spectral density of the charge excitations, lmN(q,uj), in the disordered system as the 
function of u for (a) 5 = 0.0, 0.1, 0.3, 0.5, q = 2 and (b) q= 1,2,3,8 = 0.3. 
Fig.10 The optical-conductivity is plotted for 5=0.01, and (=0.0, 0.01,0.1,0.3,0.5. 
Fig.ll The optical-conductivity is plotted for C=0.1, and 5=0.0, 0.01,0.1,0.3. 
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